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ABSTRACT 

Recent observational and theoretical studies have raised the possibility that the collimated outflows 
in gamma-ray burst (GRB) sources have two distinct components: a narrow (opening half-angle 6* j ; n ) , 
highly relativistic (initial Lorentz factor rj n > 10 2 ) outflow, from which the 7-ray emission originates, 



and a wider 



< 3#j jn ), moderately relativistic (rj w ~ 10) surrounding flow. Using a simple 



synchrotron emission model, we calculate the R-band afterglow lightcurves expected in this scenario 
and derive algebraic expressions for the flux ratios of the emission from the two jet components at the 
main transition times in the lightcurve. For viewing angles d ^ s < 9 hn we find that the contribution of 
the wide component to the optical afterglow is negligible if its kinetic energy E w is significantly smaller 
than that of the narrow component, E n , as expected for the jet core and cocoon outflow components 
in the collapsar jet-breakout model. However, if E w /E n > 1 [but the isotropic-equivalent energy ratio 
£ ; iso,w/-Eiso,n = (-B w /-Bn)(^,n/^j,w) 2 remains < 1], as expected for the decoupled neutron and proton 
components, respectively, in an initially neutron-rich, hydromagnetically accelerated jet model, then 
the narrow component only dominates the early afterglow and the wide component takes over after 
its nominal deceleration time idcc,w (typically ~ 0.1 — 1 days). Given that tdcc.w is comparable to 
the jet-break time fj e t,n of the narrow component for characteristic parameter values, the emergence 
of the wide component at idec,w may mask the jet break in the narrow component at tj C t,n, which in 
turn may lead to an overestimate of the 7-ray energy emitted by the source and hence of the required 
7-ray emission efficiency. We apply this scheme also to X-ray flash sources, which we interpret as 
GRB jets viewed at an angle 6> b s > 0j, n - Finally, we argue that a neutron-rich hydromagnetic outflow 
may naturally give rise to repeated brightening episodes in the afterglow lightcurve as observed in 
GRB 021004 and GRB 030329. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — radiation mechanisms: nonthermal 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) and their afterglows are 
commonly interpreted in terms of a relativistic outflow 
that emanates from the v icinity of a solar-mass neu - 
tron star or black hole (e.g.. lPiranl lT999: Me szaro s 2002). 
In this picture, the prompt gamma-ray emission is at- 
tributed to a highly relativistic ejecta (with an initial 
Lorentz factor 7 > 10 2 ), whereas the subsequent af- 
terglow emission in the X-ray, optical, and radio (over 
hours, days, and weeks, respectively, after the GRB) 
arises from the shock that is driven into the ambient 
medium as the ejecta sweeps up the external medium 
and decelerates. Most afterglow observations to date 
have been carried out hours to days after the GRB event, 
by which time the Lorentz factor of the afterglow shock 
has decreased to < 10. These observations have re- 
vealed the presence of achromatic breaks in the afterglow 
lightcurves of many sources, which strongly indic ate that 
GRB outflows are co llimated into narrow jets ijRhoadsl 
ll^lSarT^t^ll99l . 

Recently, however, the possibility that at least some 
GRB outflows consist of two distinct components has 
been raised in the literature. On th e observational 
side, t his possibility was first invoked bv IPedersen et abl 
(1998), who suggested that the afterglow from GRB 
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970508 could be explained in ter ms of a nar r ow je t sur- 
rounded by an isotropic outflow. iFrail et all lj2000|) sub- 
sequently proposed that the 7-rays and early (shorter- 
wavelength) afterglow emission in GRB 991216 could 
be attributed to a narrow, ultrarelativistic outflow com- 
ponent, and that the late (longer- wavelength) afterglow 
emission in this source originates in a separate wide com- 
ponent that is only mildly relativi stic. A similar picture 
was proposed for GRB 030329 by IBerger et all (|2003bl) 
and iSheth et alJ (J2003) . A two-component model was 
also suggested as an explanation of the observed re- 
brightening of the X -ray flash (XRF) source XRF 030723 
ijHuang et alJ 120041) as well as of the ap parent peak- 
energ y distribution of GRBs and XRFs ijLiang fc Pail 
,2004) and of the origin of the blueshifted optical ab- 
sorption f eatures in the spect rum of the GRB 021004 
afterglow ijStarling et alJl2005|) . 

The possibility of a two-component outflow in GRB 
sources has been independently indicated by theoreti- 
cal considerations. One can broadly divide the physi- 
cal models that give rise to such an outflow into two 
classes: models in which the separation into two compo- 
nents is an intrinsic property of the outflow, and those in 
which a narrow relativistic jet gives rise to a wider and 
slower component as it propagates through (and interacts 
with) the envelope of a progenitor massive star. One ex- 
ample of a model of the first t ype was worked out by 
iLevinson fc Eichlerl l)1993|) and Ivan Putten fc Levinsonl 
(2003): it consists of (1) a relativistic, baryon-poor jet 
driven electromagnetically along disk-anchored magnetic 
field lines that thread the horizon of a rotating black hole, 
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and (2) a subrelativistic, baryon-rich wind that is driven 
thermally from that disk. Another example is provided 
by hydromagnetically driven jets that originate from a 
neutron star or a neutron-rich ac cretion disk that form 
in the collapse of a massive star l)Vlahakis et alJ 12003'!. 
In this case the neutrons decouple at a moderate Lorentz 
factor while the protons continue to be accelerated and 
collimated by the electromagnetic forces, giving rise to 
a narrow, highly relativistic proton component and a 
wider and slower neutron component (which, after decou- 
pling, is transformed into a moderately relativistic proton 
component through neutron decay). Examples of mod- 
els of the second type include jet-induced core-collapse 
supcrnovae, wherein collimated high-velocity jets cause 
the envelope of the progenitor massive star through 
which they propagate to be ejec ted with subrelativisti c 
speeds and an oblate geometry (Khokhl ov et al.1 11999'1. 
and the collapsar model, in which the outflow result- 
ing from the jet breakout through the progenitor star's 
envelope is predicted to consist of a highly relativistic 
jet c ore and a moderately relativis tic surrounding co- 
coon l)W. Zhang et alJl200l see alsoE amirez-Ruiz et aLl 
2002). 

In this paper we focus on two-component GRB out- 
flows in which both components initially move with rel- 
ativistic speeds and therefore end up contributing to the 
optical afterglow emission. Accordingly, we adopt as rep- 
resentative examples the hydromagn etically accel e rated , 
initially neutron-rich jet model of iVlaha kis^^d] l|2QQ3j) 
and th e collapsar jet-breakout model of IW. Zhang^taLl 
(2004). According to the numerical simulations of the 
latter authors, the narrow component in the collapsar 
model has a Lorentz factor i] D > 100 and an opening half- 
angle 6>j jn ~ 3° — 5°, whereas the corresponding quantities 
for the wide component are ry w ~ 15 and 9^ w ~ 10°, re- 
spectively. The characteristic Lorentz factors in this sce- 
nario are very similar to those {rj n ~ 200 and r/ w ~ 15) 
in the representative ne utron-rich hydromagnetic model 
of lVlahakis et alJ (J2003). 3 However, in contrast with the 
collapsar model, in which the highly relativistic jet com- 
ponent is in general more energetic (typically by about 
an order of magnitude) than the cocoon material, the 
asymptotic kinetic energy i? w of the wide component in 
the hydromagnetic jet model typically exceeds the cor- 
responding energy of th e narrow component (E w ss 2 E n 
in the fiducial model of Vlahakis et al. 2003). 

Our goal in the present work is to examine some of 
the observational properties of two-component GRB out- 
flows. In particular, we calculate (§ 01 the approximate 
optical afterglow lightcurves that are produced by the 
shocks that the two jet components would drive into the 
ambient medium. We then argue (§ EJl that outflows 
of this type may have significant general implications to 
our understanding of GRB and XRF sources. Our con- 
clusions are given in § 0] 

2. MODEL AFTERGLOW LIGHTCURVES 

A simple jet structure is assumed in this work, con- 
sisting of a narrow and initially faster component and 
a wide and initially slower component. Each compo- 

3 In the simplified model used by [Vlahakis ct al. 1 2003}), the 
value of 0j jW could not be calculated exactly; in this paper we 
assume that it can be as large as ~ 3 8j, n - 



nent is assumed to be uniform within some finite open- 
ing angle and to have sharp edges. This two-component 
model should not be regarded as a limiting case of the 
structured " universal" jet models discussed in the lit- 
erature (e.g-lRossi et al.ll2002t iZhang fc Meszarosll2002l 
IB. Zhang et al .1120041) . In the latter models, all jets are 
nearly identical and their injected energy per unit solid 
angle has a power-law or a Gaussian dependence on the 
polar angle 9 (measured with respect to the jet axis). 
In contrast, in the scenario that we consider the wide 
component does not contribute to the 7-ray emission 
and hence cannot be a part of a traditional universal-jet 
model. In the uniform, sharp-edged jet picture the open- 
ing angles of the two outflow components are not invari- 
ant from source to source although their values may well 
be correlated. It is, however, conceivable that each of the 
two components is structured and that the jet is "uni- 
versal" in the sense that this structure varies little from 
source to source. The basic implications of a structured 
two-component jet model would be similar to the ones 
that we discuss in § [3] in the context of a sharp-edged, 
uniform outflow, although some of the details may be 
different and would depend on the specifics of the angu- 
lar distribution of i? an d n in each component (s ee, e.g., 
Kumar fc G7anotjl2003l and IB. Zhang et"aTll2004^ . 

In our simple treatment the interaction between the 
two jet components is neglected. This assumption can 
be justified even after the narrow component's jet-break 
time (see eq. below) , when the effects of side ways ex- 
pansi on could in principle become relevant (e.g.. lRhoadsl 
1999), in view of indications from recent numerical sim- 
ulatio ns (e.g., iKumar fc Granotl 120031: iCannizzo et alJ 
2004) that in practice there is relatively little lateral 
spreading so long as the jet is at least moderately rel- 
ativistic. To further simplify the discussion we only con- 
sider the case of a uniform external medium (of number 
density n = hq cm~ 3 ) and we neglect the possible ef- 
fects of radiative losses on the hydrodynamic evolution 
(which may affect the early afterglow during the period 
of fast cooling if the fraction e e of the internal energy in 
electrons behind the afterglow shock is not <C 1). The 
narrow and fast jet component has an initial Lorentz 
factor T] n (> 10 2 ), a half-opening angle (9j in , and a ki- 
netic energy (at the beginning of the afterglow phase) 
E n , while the wide and slow jet component is character- 
ized by 7/ w (~ 10), #j, w (> 0j )n ), and E w . In the fol- 
lowing, the subscripts 'n' and 'w' will denote the narrow 
and wide jet components, respectively. The ratio of the 
true energy E and the isotropic-equivalent energy E- mo 
is given by the beaming factor f b = l — cosflj « 9?/2. 

Thus, i?iso,w/i?iso,n = (#j, n / '#j,w) 2 E w / ' E n . 

The emission from each outflow component is calcu- 
lated separately. For the early afterglow (while the 
reverse shock is still present) we use the results of 
ISari fc Piranl (|1999albf) . For t he emission during the sub- 
sequent se lf-similar evolutio n ijBlandford fc McKeell976() 
we follow ISari et alJ l|1998fl . an d for the post-je t-break 
emission we use the results of ISari et all ^99). The 
typical synchrotron frequency v m , the cooling frequency 
is c , and the peak flux -Fi, jma x, of the emission from the 
shocked external medium behind the forward shock are 
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given by 



v m = l.l x 10 1 % 2 e2_ 1 e^ 2 _ 1 nJ /2 ( 7 /300) 4 Hz , (1) 

C 2 (7/300)- 4 Hz , (2) 



1/2 3/2 n-2 + 3 



220 £b.-i«o ^i:,28^(7/300) B M Jy 



(3) 



l|Sari fc Piranlll999b|) . where g = 3(p - 2)/(p - 1), p 
is the power-law index of electron energy distribution 
(dN e /dje oc 7~ p ), f = f s sec is the observed time, 7 is 
the Lorentz factor of the shocked fluid, es is the fraction 
of the internal energy behind the shock in the magnetic 
field, Dl is the the luminosity distance, and Qi = Q/(10 l 
times the c.g.s. units of Q). 

The interaction of the jet with the ambient medium ini- 
tially drives a reverse shock into the GRB ejecta, which 
decelerates the ejecta. When the reverse shock is New- 
tonian, or at most mildly relativistic, then 7 » r\ over 
its entire duration and the energy given to the swept-up 
external medium (of rest mass M) is j 2 Mc 2 ~ i] 2 Mc 2 . 
As this energy approaches E, the original kinetic energy 
of the ejecta, after 



fde 



Rdcc 

2cr] 2 



= 0.49 



E; 



iso,52 

n 



1/3 



-8/3 



days , (4) 



significant deceleration must occur. For f > fdcc most 
of the energy is in the shocked ext ernal medium and a 
self-s imilar evolution is established ( Blandfor d fc McKed 
Since i] w is assumed to be rather small (~ 10), 
idoc,w (~ 0.5 days) is much larger than the duration of 
the GRB. The refore, the ej ecta is alw ays in the "thin 
shell" regime ISari fc Pira.nl 179951 Islrll[l99l, In this 
case the reverse shock is initially Newtonian. It is nat- 
ural to expect some variation in the initial Lorentz fac- 
tor, Ari ~ T]. For a thin shell, this causes the ejecta 
shell to start spreading long before the reverse shock fin- 
ishes crossing the shell, which in turn causes the reverse 
shock to become mildly relativistic before the crossing 
has ended. In this paper we concentrate on the opti- 
cal emission from the forward shock, which in the case of 
the fast component would dominate the optical flux from 
the reverse shock after ~ 1 3 sec for typical parameters 
ijKobavashi fc Z hang 2003). In the case of the slow com- 
ponent, the contribution from the forward shock should 
typically dominate the optical flux at all times, with the 
emission from the reverse shock maki ng a significant con - 
tribution only in the radio band fe.g.. lPiran et alJl2004fl . 

When 7 drops to ~ O^ 1 the edge of the jet becomes 
visible and sideways expansion may become noticeable. 
These two effects cause a break in the lightcurve at 



ijot = 0.25 £7: 



1/3 



-1/3.8/3 



iso,52"0 



0> days , 



(5) 



with the former effect evidently responsible for most of 
the steepening if the Lorentz factor is not too close to 1. 
Expressed in terms of the true energy, the jet-break time 
is 

f jet = 0.66 El{\- 1/3 2 days . (6) 



In the early afterglow, at t < fdcc, 7 ~ V l|Sari fc PirarJ 

Il999albl). At W < t < fjet, 7 ~ 0f HVtjet) -878 

ijSari et al.lll998(l . At t > fj et (and before the nonrela- 
tivistic transition time ^nr), we have 7 ~ O^ 1 (t/t^t)^ 1 / 2 



assuming rapid lateral expansion ijSari et al.1 [l999). 
Therefore, the temporal scalings of the break frequen- 
cies and peak flux arc given by 



f , f *C fdcc , 
J' ;l , CX7 4 OC { f _3/2 ,f d cc < t < fjet , 
t > fjet j 

t , f < fdcc , 
/'< X7~ 4 f~ 2 CX { f _1 / 2 ,fdcc < f < f 
f t > fj t 



t 



-2 



jet , 



f , t <C t dec , 
4-0 



OC7 f OC < f U , fdcc < f < fjet , 
t , t > fjet • 



(7) 



(8) 



(9) 



In the limit of negligible sideways expansion after the jet- 
break time, the time dependence of 7 does not change at 
fjet, so the behavior described by the second line in equa- 
tions and (JSJ) continues to hold also for t > fj c t- How- 
ever, the maximum flux in this case is still reduced by a 
factor (#j7) 2 (representing the ratio of the jet area to the 
beaming cone area) as f increases above fj e t, so f _1 in the 
last line of equation © is replaced by f _3 / 4 . For practical 
applications the qualitative behavior of the lightcurve in 
this limit is very similar to that in the limit of rapid lat- 
eral expansion, and since the expressions given in equa- 
tions 0-© are the ones commonly used in the litera- 
ture, we continue to employ them in this work. 

The transition time f from fast cooling to slow cooling 
(when v m — v c ), and the times f m and t c when v m and 
u c , respe ctively, pass by the observed frequency v, are 
given by (|Sari et a,l. 111991 

f = 0.55 g 2 E 1S0i52 e 2 -i n o hr , (10) 



f m =o.36. 9 4 / 3 4/ 3 52 ,r 5 2/3 ^ 4/3 - 1 1 

n-1 ..-2,-3 -2 



C e,-1 C B 



(11) 

f c = 0.17 ^vrsV-i^ 3 hr • ( 12 ) 

The transition frequency, vq, defined by f m (fo) = 
v c (to), is given by 

,-5/2 -3/2 



5.5 x 10 14 g-'Er^e- 1 ^ 



Hz 



(13) 



J iso,52 c e,-1 C B,-1"-0 

These transition times together with fd cc and fj e t sep- 
arate the time domain into several segments. At each 
time segment, the flux is derived by comparing the ob- 
served frequency to v m and v c to determine the appro- 
priate spectral behavior; one also makes sure to use the 
correct dynamical behavior for the given time segment. 
The explicit expressions for the flux contributions of the 
two outflow components at the different time segments 
are presented in Appendix lAl 

The ratio of the deceleration times of the two jet com- 
ponents is 

1/3 / \ -8/3 

fe) ' <»» 

10 3 . The deceleration time of 



fdc. 
fde 



which for rj n /r] w ~ 10 is 
the slow component is much larger than that of the fast 
component, so a bump would show up in the decaying 
lightcurve of the fast component due to the emission of 
the slow component if -F„ iW > F v n at f = f dcc.w ■ 

The flux ratio (i 7 i/, w /^,n)t=t d cc w depends on whether 
the slow/wide component decelerates before the jet- 
break time of the fast/narrow component or not (i.e., 
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, , w and £jet,n), since 7 has 



on the relative ordering of id- 
different time evolution indices before and after jet break 
From equations (@J and © , the ratio of these two times 
is given by 



A 1 



Eu 



1/3 



-8/3 



^jct,n 

where A x = A 2 C^ t 3 and A 2 = C t (t > t dcc )/C t (t dcc ) ~ 
(with the coefficients C ]C t and Ct defined from 7(ij t) 
'jet/^j and t = R/Ct"f 2 c). In this work we assume Cj 



(15) 



C 



jet 



1, Ct(t dcc ) = 2, and C t (t > t dcc ) = 4, for which A\ 
A2 = 2. However, in the figures of this paper we calculate 
the flux by using Ct — 4 at all times to ensure continuity 
of the plotted lightcurves; our algebraic expressions for 
the component flux ratios may therefore yield values that 
differ somewhat (by a factor < 2) from those implied by 
the presented figures. 

From equations (0J, and the scaling of 7 with time, 
the peak-flux ratio of the two components at t deCiW is 

-Eiso 



F„ 



(^7) ( wr) tdcc ' v 



< t. 
> t 



jct,n 



(16) 



jet.n 



The flux ratio under consideration also depends on the 
frequency range within which the observed frequency v 
is located at idec,w for each of the two outflow compo- 
nents. For typical parameter values, idcc,w > io,n> so 
the narrow/fast component is in the slow cooling regime, 
^c,n(idec,w) > ^m,n (idoc,w) • We also expec t the wide/slow 
component to be slowly cooling (see eq. IA5| ). 

Based on the discussion in Appendix ^ the R-band 
observation frequency ^ = 5x 10 14 Hz exceeds both 



,(td 



and Vm, w (t dec ^ n ) for typical parameter val- 



ues. However, can be either larger or smaller than v c 
for both components. There are thus four relevant cases: 
If v m (tdcc. w) <v< v c (t dcc w ) for both components, then 

(p-i)/2 

= h 



F, 



v,max,w 



F 



i/.max.n 



.4 



3(p-l)/4 



(fe 



(p+3)/4 



^dec,w ^jet,n 



A 



(zH 



3)/4 j4 3(p-l)/4 



X (7] u 



2(p+3)/3 td 



(p+3)/3 



(17) 



> t 



jet, 11 



This parameter regime applies to all the numerical ex- 
amples presented in § |3 If v > v c (t<icc,w) > ^m^dcc.w) 
for both components, then 



F„ 



F 



v,max,w 



F 



1/2 



A 



3(p-2)/4 



(fe) 



(p+2)/4 



(P"l)/2 



tde 



= h 



£jet,n j 

(18) 



A 



OH 



2)/4 j4 3(p-2)/4 



(fc) 

2(p+2)/3 f(J 



If V > V c , n (t d 

cc,w ^ j 



1 

X (»7«)0j,n 
^m,w(^dcc,w) <C ^c,w 

1/2 



(p+2)/3 



= A — 



= / 2 ( 



(*dec,w), then 
1/2 



(19) 



If fm,n(*dec,w) < ^ < ^c,n(idcc,w), V > ^c,w(^dcc,w), then 



= /l 



-1/2 



-1/2 



(20) 



We find that, at ij e t,w, Vr > and v c > v m for both 
components over most of the characteristic parameter 
ranges (see AppendixEJ. Given also that the flux evolu- 
tion after the jet break is the same for both components, 
it follows that the flux ratio after the two components 
have undergone a jet break is the same as that at ij e t,w 
From equations ©, ©, and 0, the peak- flux ratio of 
the two components at t > tj C t,w is 4 



F 



t>t u 



-^iso,w 



^jet,w 



jet,n 



E ^ 4/3 

(21) 

Therefore (specializing to the case tdcc.w > ijet,n) 5 the 
flux ratio at tj e t, w assumes the same form as the flux ra- 
tio at £dec,w after rj w is replaced by l/0j C t, w and all the 
characteristic frequencies are evaluated at tj e t,w Denot- 
ing by fi the same flux ratios as fi for i = 1, 2, but 
evaluated at ij c t,w instead of at t dcCjW , we obtain the fol- 
lowing simple results: 

1 _ 



p 



E w \ (p+3)/3 



and 



± i/,w 



~E n 



(p+2)/3 



(22) 



(23) 



Another interesting quantity is the ratio of the fluxes 
from the two outflow components at the corresponding 
jet-break times, / = i^, w (tj e t, w )/^/,n(*jet,ti)- We find 



fi 



^w,w (^jet,w) 



-^y,n(^jet,n) 



-2p 



(24) 



and 



F„,n(t 



2/3 



2p 

^] .(25) 



So far the calculations correspond to on-axis observers, 
i.e., # bs = 0. The results, however, apply to the 
entire range of observation angles 6* b s ^ 1.5 0j, as 
demonstrated with the help of more realistic jet models 
(|Granot et al.ll200^) . For off- axis observers with viewing 
angle # b s ^ 1-5 #j the afterglow emission peaks at t$, 
the time when 7 



'obs : 



D 



'obs 



'■jet 



(26) 



In the simple jet model that we use IRhoadslll99flt ISari et alJ 
119991) . the jet dy namics becomes ind ependent of its initial opening 
angle at t > tj ct I Granot ct al. 200.3). Therefore, the lightcurves at 
any given viewing angle 6 ^ s depend only on the true energy E of 
the outflow (in addition to the ambient density and microphysical 
parameters e e and eb> which are assumed to be the same for the 
two jet components). Thus, the flux ratios at a given observed time 
and frequency after both jet breaks (and both deceleration times) 
depend only on the ratio of their true energies, E„/E n . 
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ijNakar et alJ200^1 . where the model-dependent factor B 
is of the order of unity (we adopt B = 1 for the numerical 
estimates in this work). The lightcurve for t > tg is 
similar to the on-axis lightcurve. The maximum flux at 
tg strongly decreases with increasing viewing angle, as 
Fna^ ujOobs) oc 0~£* (see eqs. [9] and [10] of iNakar et~aTl 
2002). As we discuss below, for 1.5dj >n < o bs ^ 1.5 0j, w 
it is possible for the wide component to be visible first 
and for the narrow component to become dominant (at 
least temporarily) later. The estimate of equation 115(1 
is now replaced by 



^dec ' 



B \ E n J 9 bs 



'j,w^7w j 



-8/3 



(27) 



The emission from the wide and the narrow outflow com- 
ponents peaks at idcc,w and tg >TL , respectively. In evalu- 
ating the effect of these peaks on the overall lightcurve, 
it is useful to consider two flux ratios: the quotient 
f a of the first and second peak components, given by 
F v , w (td e c,w)/F v ,n(t0,n) if tg,n > <dcc,w and by its in- 
verse if ig, n < idcc.w, and the ratio f b of the pri- 
mary and secondary flux contributions at the time t = 
max{£dcc,w, te, n } of the second peak. The ratio f a deter- 
mines which component dominates the overall lightcurve, 
whereas f b indicates whether the secondary component 
can play a role in the late afterglow. As is the case at 
tjct.w, typically the R-band observation frequency at tg >n 
exceeds u m and u c > v m for both outflow components. 

If idec.w > te,n then f b (which in this case is the flux 
ratio of the wide component to the narrow component at 
^dcc.w) is the same as the ratio / obtained in the on-axis 
case for idcc,w > £jet,n- The difference from the on-axis 
case is that the flux of the narrow component peaks at a 
later time (ie, n ) and is smaller. In this case 

f a = fb (^n_Y (28) 



or, written down explicitly. 
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for the cases where v < v c and v > v c , respectively. In 
order for a bump associated with the wide component to 
become visible during the late afterglow it is necessary 
for f b to exceed 1. This condition requires E W /E D > 

A -3/4 A -9( P -2)/i( P+ 2) M ^ 2 (gee ^ jjjj jjg^ 

If f a is also > 1 then the wide component dominates the 
entire lightcurve, with the narrow component possibly 
becoming visible as a bump on the curve's rising branch. 
If, however, both f a and f b are < 1 then the narrow 
component would dominate at all times. 
When i doc , w < tg in < £jet,w, f b is given by 
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whereas f a is given by the inverse of the expression (eq. 

for idcc.w > tg : n- In this case the condition f b > 1, 
which requires E n /E w > -B(#obs/#j,w) 2 (see eqs. and 
|32|). corresponds to the narrow component dominating 
the late afterglow. If also f a > 1, then the narrow com- 
ponent dominates the entire lightcurve, with the wide 
component possibly becoming visible as a bump during 
the curve's initial rise. Conversely, when both f a and 
f b are < 1, then only the wide component's afterglow 
emission would be visible. 

If the observer is located outside the solid angle sub- 
tended by the wide outflow component (i.e., # b s ^ 
1.5 ^j,w) , then the flux contributions from the wide 
and the narrow components peak at i# iW and tg tIl , re- 
spectively. The ratio of these times is tg lW /tg :a = 
(Ew/En) 1 / 3 , which is independent of 9 b s - Thus, when 
E w > E D then tg yW > tg^ n , /" = E W /E D , f% = 

(E w /E n ) 2 / 3 , and f b = f. However, if E w < E n then 
to.w < t e , n , fi = E n /E Wl ft = {E n /E w ) 2 ' 3 , and 
f b = All of these cases have f a > 1 and f b > 1, 

which means that the more energetic outflow component 
dominates the lightcurve for t > max{te iW , tg^ n }. 

3. APPLICATIONS 

Our results have potentially significant implications to 
the interpretation of GRBs and XRFs. We consider 
these two applications separately, although our discus- 
sion makes it clear that they could be related under a 
unified picture of these sources. 

3.1. GRB Afterglows and Source Energetics 

We choose as fiducial parameters rj n — 200, ?/ w = 15, 
0j, n = 0.05 rad, 0j, w /0j,n = 3, andp = 2.2. For this choice 
the ratio idcc.wAjct.n of the deceleration time of the wide 
component to the jet-break time of the narrow compo- 
nent (eq. [T5|l ranges between 3.0 and 1.4 as E^/E^ de- 
creases from 3 to 1/3. The near-coincidence of these two 
time scales can have interesting observational ramifica- 
tions. In particular, if the flux ratio {F v ^ / 'i^, n )t=t dec w is 
close to 1, then the presence of a break in the narrow com- 
ponent at ijot.n may be masked by the rise in the flux from 
the wide component that occurs as idec,w is approached. 
For the adopted parameters, one finds from equations 
l(T7 |) -lJ2U jl that, in fact, (F„, w /F v 

,n)i=t d cc,w ^ 1 SO long 

as E W /E D is > 2, but that the flux ratio becomes -C 1 
for low values of the wide-to-narrow injected energy ra- 
tio. For example, in the case described by equation (|17(l . 
(F UtVr /F v , n )t= tdoc iW = 1.9, 0.3, and 0.04 for E w /E n = 3, 
1, and 1/3, respectively. For t > £dcc,w > tjet,n the 
flux from the wide component decreases with time no 
faster than t - ( 3 P~ 2 )/ 4 (= i -115 ) or even increases with t 



(for i d < i < min{£ 



see eqs. |A"fj ] -|ATT ] ) 



m,w j ^jet,w 

whereas the flux from the narr ow component decreases 
steeply (as £~ p ; eqs. EHUD' Th us, even if F v w is 
still < F v ,n at £ = i doc ,w (i.e., /i < 1), when E w > E n 
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Time [days] 

Fig. 1. — R-band afterglow lightcurve from a two-component 
jet. The contribution of the narrow component, wide component, 
and their sum is represented by the dashed, dash-dotted, and solid 
curves, respectively. The total outflow energy is assumed to be 
constant, E v + E n = 10 51 ergs. From top to bottom, the 3 panels 
correspond to E w /_E n = 3, 1, and 1/3, respectively. The other 
parameters are the same for all panels: rj n = 200, n w = 15, 6\ v = 
0.15, 6»j n = 0.05, n = 1, e„ = 0.1, t B = 0.01, p = 2.2, 'and 
£>L,28 = 1. 



it will become the dominant contributor to the total af- 
terglow flux soon thereafter. (In the case described by 

eq. [TZ|, this will occur for i/tdcc,w > f\ 4 ^ p+2 \) As 
illustrated by the model lightcurve plotted in Figure QJi, 
a clear signature of a jet break in the lightcurve would 
occur, under these circumstances, only at t = ij c t,w- 

The possibility that the jet break in the narrow outflow 
component remains unobservable could have important 
consequences for the infer red energetics of GRBs. Recall 
that it has been found ijFrail et alJ 120011 iBloom et al.l 
|2003|1 that, when the isotropic-equivalent 7-ray ener- 
gies of a sample of GRBs are converted into true en- 
ergies by using the beaming factor inferred from the ob- 
served jet break, then the resulting values cluster nar- 
rowly about E 7 « 10 51 ergs. On the other hand, the 
isotropic-equivalent kinetic energies of GRB outflows, as 
estimated from dynamical and spectral modeling of the 
associated afterglows, yield — after being corrected by 
the jet break-inferred beaming factor — a narrow dis- 
tribution of true kinetic energies (for 2 oppositely di- 
rected jets) at the beginning of the afterglow phase that 
is center ed onE k 10 51 er gs ()Panaitescu fc Kumarl2002T: 
see also lYost et al.ll2~003|) . A similar result is obtained 
when the afterglow X-ray luminosity is used as a sur- 
rogat e for the isotropic- equivalent outflow kinetic en- 
ergy ijBerger et al.ll2003ajL The X-ray emission typically 
peaks during the early afterglow, so the kinetic energy 
estimated on the basis of the X-ray luminosity (and con- 



ventionally evaluated at t — 10 hr) likely corresponds 
to that of the narrow outflow component. However, the 
X-ray-based isotropic-equivalent kinetic energy is typi- 
cally found to be smaller than the isotropic-equivalent 
kinetic energy inferred from the spectral and dynamical 
modeling of the overall afterglow (which, for E w > E n , is 
dominated by the contribution of the wide component). 
In the context of the two-component jet model this im- 
plies that the wide component should dominate also at 
early times, which is clearly inconsistent. The X-ray- 
based isotropic-equivalent kinetic energy typically also 
turns out to be smaller than the isotropic-equivalent 7- 
ray energy, which in a two-component model that asso- 
ciates the 7-rays with the narrow core would be difficult 
to reconcile with the internal-shock scenario of GRBs (see 
discussion in the next paragraph). It is therefore quite 
possible that the X-ray-based deduction systematically 
underestimates the true kinetic energy in the narrow out- 
flow component. 

The approximate equality of the inferred values of 
Ey and E has been given several different explanations; 
here we focus on its interpretation i n the context of 
the in ternal-shock scenario for GRBs i|Rees fc Meszarosl 
1994), which has been successful at accounting for 
the observed variab ility properties of the bursts (e.g., 
iNakar & Piranll2002|) . In this picture, the 7-ray emis- 
sion originates in shocks that form in the collisions of 
"shells" that are injected with vari able energy and /or 
mass at the origin. It was shown ( Be loborodovl 1200(1 
iKobavashi fc Sarill200lD that E~ ( /E can in principle be 
~ 1 in this case if the following conditions are satis- 
fied: (1) the spread between the minimum and maxi- 
mum initial Lorentz factors of the shells is large enough 
(7i,max/7i,min S 10); (2) the distribution of initial 
Lorentz factors is sufficiently nonuniform (one possibil- 
ity being that log7i, rather than 7i, is distributed uni- 
formly); (3) the shells are approximately of equal mass 
and their number is large enough (> 30), and (4) the frac- 
tion of the dissipated energy that is deposited in electrons 
and then radiated away is sufficiently high (e e > 0.3), 
with a similar constraint applying to the fraction of the 
radiated energy that is emitted as gamma-rays. If any 
of these conditions were violated to a significant extent 
then the implied magnitude of E 1 /E could decrease to a 
value well below 1. 

We do not at present have independent information 
about the nature of GRB outflows to verify that the 
above conditions are indeed satisfied in the majority of 
sources, as would be required for consistency between 
the internal-shock model and the inferred distributions 
of and E. It is, however, worth noting that these con- 
straints can in principle be alleviated if the outflows cor- 
respond to two-component jets with E w > E n . 5 This is 
because the inference E 1 jE ~ 1 is based on the assump- 
tion that the solid angle used to convert the isotropic- 
equivalent energy into a true energy is the same for the 
7-rays and the afterglow radiation. If E w > E n then 
most of the afterglow radiation will be emitted from the 
wide component and the value of E will be appropriately 
inferred from Ei so using the opening half-angle of the 
wide outflow component, 0j,w- However, as the prompt 

5 This possibility was originally note d in a talk at the 2003 GRB 
meeting in Santa Fe; seefKonigl (2004|). 
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high-energy emission in this picture originates in the nar- 
row outflow component, the conversion of the isotropic- 
equivalent 7-ray energy into E 7 must be done using the 
opening half-angle of the narrow component, 0j. n . If, 
as discussed above, the jet break in the narrow compo- 
nent is not observationally discernible and the outflow 
is mistakenly interpreted as a single-component jet with 
an opening half- angle 0j, w > then E 7 will be overestimated 
by a factor ~ (^j,w/^j,n) 2 (= 9 for the fiducial values 
adopted in this paper). The actual magnitude of £L in 
this case could thus be well below the value inferred on 
t he basis of a sing l e-comp onent jet model. 

iKumar fc Piranl l)2000bt) proposed that the ejected ma- 
terial in GRB outflows exhibits strong angular fluctua- 
tions, forming "patchy" shells. In this case the conversion 
from an isotropic-equivalent to a true energy also involves 
a smaller effective solid angle for the 7-ray emission than 
for the afterglow radiation. This situation could in prin- 
ciple be distinguished from the two-component outflow 
scenario discussed in this paper through some of the 
spec ific predictions o f each o f these models. For exam- 
ple, IKumar Ez Piranl l)2000b(l argue that a patchy-shell 
outflow could exhibit large temporal fluctuations (with 
a progressively decreasing amplitude) during the early 
(minutes to hours) afterglow, whereas the results derived 
in this paper indicate that the afterglow lightcurve pro- 
duced by a two-component jet with E w > E n might tem- 
porarily depart from a simple power-law decay around 
t ~ £dcc,w (over a timescale of hours). 6 It is, however, 
conceivable that the ejected shell may be patchy even if 
the outflow has more than one component. 

As we have shown, the possibility that E 1 is overesti- 
mated in a two-component outflow can only be realized if 
E w exceeds E n . However, in order for the high-efficiency 
requirement on the emission from internal shocks (and 
the corresponding conditions listed above) to be relaxed, 
the ratio E w /E n cannot be much greater than 1. Specifi- 
cally, the kinetic-to-radiative energy conversion efficiency 
of the narrow outflow component, £ n = E 1 j(E 1 + E n ), 
is determined by the ratio Ey/E n , which is overesti- 
mated by the factor (E n /E w )(8 jtW /8^ n ) 2 = E is0tn /E is0tW 
if E w E w . 7 Thus, to have any reduction in the required 
efficiency, i?i so ,n/-E ; iso,w must exceed I. The above two 
constraints can be expressed as a double inequality on 
the ratio of the component kinetic energies: 

1 < E w /E n < (^w/^n) 2 ■ (33) 

The condition E w /E n > I implies that the wide compo- 
nent dominates the afterglow emission at late times (see 
eqs. |21 | -|23 | ). whereas the requirement £^ S o,n/-^iso,w > 1 
implies that the narrow component dominates at early 
times (see eqs. [T?)]-|271]). 

A two-component outflow with E w /E n > 1 arises 
naturally in the (initially) neutron-rich, hydromagnet- 
ically accelerated jet scenario (see § QJ. In contrast, 
in the two-component outflow investigated in the con- 
text of the collapsar model, the ratio E w /E n is typ- 
ically <C 1 and hence (see Fig. ^) the optical after- 

6 As discussed below, refreshed shocks are another likely source 
of lightcurve variability (on timescales of hours to days). 

7 If the kinetic energy inferred from the afterglow fitting in fact 
corresponds to the sum of the contributions from the wide and 
narrow components, then the factor E n /E w in the overestimation 
expression is replaced by 1/(1 + E w /E n ). 




10 * 10 10 u 10 10^ 

Time [days] 

Fig. 2. — Similar to Fig. except that _E W = 2 E n and 0j jW = 
2 0j, n - The top and bottom panels correspond to 0j n = 0.05 and 
0.1, respectively. The other parameters are the same as in Fig. HI 



glow emission from the wide (cocoon) component would 
generally remain undetectable at all times. For the 
adopted fiducial parameters and assuming E w /E n =0.1, 
£dcc,wAjet,n = 0.96 and, for the case described by equa- 
tion (F VjW /F Vin ) t=tdm w = 5 x lfT 3 . Note, how- 
ever, that the cocoon afterglow emission might dominate 
at early times at submillimcter wavelengths and that 
there may also be a signature in the early optical after- 
glow of the collision between the expanding cocoon and 
the decelerating head of t he narrow outflow component 
(jRamirez Ruiz et al.l l2002) . 

The plots in Figure demonstrate that the optical 
afterglow lightcurve from a two-component jet departs 
from a simple power-law decay when the flux contri- 
butions from the two components become comparable, 
which for a jet with E w > E n typically occurs around 
the deceleration time of the wide component. As the 
wide component gradually takes over from the narrow 
component to become the dominant contributor to the 
flux, the lightcurve exhibits a concave "flattening" (if 
tdec.w > tjet,n> see Fig. \Bp) or a convex "bump" (if 
tdec.w < f jot,n; see Fig. 0>) of duration At ~ t. The pres- 
ence of this feature may be hard to discern in practice be- 
cause of insufficiently dense time coverage or the interfer- 
ence of other factors (the emission from the reverse shock, 
ambient density inhomogeneities, refreshed shocks, etc.) 
that could cause a similar behavior. Nonetheless, there 
are already several potential candidates for this feature 
among observed afterglows. For example, GRB 970508 
exhibited a pronounced brightening between ~ 1 and 
~ 2 days, after w hich it followed an ap proximate power- 
law decay (e.g., IPedersen et al.lll998j) . However, be- 
tween ~ 0.1 and ~ 1 day the lightcurve was constant or 
slightly declining with time, a behavior that is not repro- 
duced by our simple two-component jet model. 8 A ma- 

8 In our picture, the brightening would be caused by a wide 
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jor brightening event wa s also recorded i n GRB 021004 
around ~ 0.1 days (e.g.. iFox et al]l2~003D . In this case, 
the lightcurve assumed a flat form between ~ 0.02 and 
~ 0.1 days and a power-law decay ind ex of ~ 1.2 immedi- 
ately thereafter ljUemura et all2003fl . which is consistent 
with the behavior expected in the two-component model 
(see Figs. ^ and|2t; if this brightening is indeed associ- 
ated with the emergence of the wide-component emi ssion 
then the inferred power-law index implies, by eq. |A6| . 
p w 2.2 for this component). The identification of the 
brightening time with idcc,w suggests that the wide com- 
ponent in this source had a comparatively high value of 
r] w (see eq. g]). 

An alternative interpretat ion of the early brightenin g 
in GRB 021004 was given bv lKobavashi fc Zhanel {2003), 
who attributed it to the emission of the forward shock 
taking over from that of the reverse shock. It is, how- 
ever, worth noting that GRB 021004 exhibited a second, 
less pronounced brightening at t w 1 day and possibly a 
third one at t « 3 days, and that it has been suggested 
that all these events may have a similar physical ori- 
gin — either a variable external density ijLazzati et al.l 
.2002) or energy fluctuations that, in turn, could arise 
either from variable injection at the source (refreshed 
shock s) or from a patc hy angular structure of the out- 
flow ijNakar et al]l2003l) . Interestingly, a refreshed-shock 
scenario is a natural feature of the hydromagnetic. ini- 
tially neutron-rich jet model of lVlahakis et alll|2003|) . In 
this picture, the decoupled neutrons that constitute the 
wide outflow component decay into protons on a distance 
scale > Rp = 4 x 10 14 (7y w /15) cm, which is likely larger 
than the scale over which many of the shell collisions in- 
voked in the internal shock model for GRBs take place. 
Shell collisions may well give rise to the 7-ray emission 
from the narrow (proton) outflow component, but since 
they cannot take place inside the wide component be- 
fore the neutrons are converted into protons, all the fast 
shells that overtake slower shells at R < Rp would be- 
come arranged in a sequence wherein the faster shells are 
in the front and the slower ones lag behind and remain 
closer to the origin. Furthermore, the radius at which 
the neutron shells are arranged in this way is smaller by 
a factor of ~ (i] n /i] w ) 2 compared to the radius of the 
internal shocks, and is thus <C Rp. The neutron shells 
can therefore pass through each other with very little 
interaction between them while ordering themselves ac- 
cording to their velocities. After the decay into protons 
the wide outflow component would start sweeping up the 
ambient mass, which would cause it to decelerate. Under 
these circumstances, the slower shells that had been left 
behind would overtake the decelerated front shell, lead- 
ing to a pileup as progressively lower-n shells arrive at 
correspondingly later times. The wide-component after- 
glow would then assume the form of a repeatedly reen- 
ergized shock, with the energy injection occurring quasi- 
continuously at first and then possibly tapering off as the 
slowest shells finally arrive at the front-shock location. 
This picture is broadly compatible with the observations 
of GRB 021004: the large initial brightening may be in- 
component with E w > E n that is observed at o h s 5 -, $j, n an( i peaks 
at t = t^cc w In contrast, Panaitcscu ct a!,] 119981) attributed the 
flux rise to a narrow jet seen outside its opening half-angle and 
suggested that the flux at earlier times could be produced by a 
wide (essentially isotropic) component of lower energy. 



terpreted as the quasi-continuously energized early after- 
glow emission from the wide outflow component, and the 
subsequent rebrightenings may be attr ibuted either to 
collisio ns with late-arriving shells (e.g., [Kumar fc Piranl 
2000a) or to a patchy angular structure (for which other 
aspects of the a fterglow p rovide independent s upport; see 
INakar fc Piranll2003l and iNakar fc Orenll2004[) . 

Similarly to GRB 021004, GRB 030329 manifested 
a significant brightening in its optical lightcurve (at 
i ~ 1.5 days), followed by several less pronounced re- 
brighte nings (at t ~ 2.6, 3.3, and 5.3 days, res pec- 
tively). iGranot et all ((2003) and lPiran et~aTl l|2004D ar- 
gued that these brightening episodes can be interpreted 
in terms of refreshed shocks in a single-component jet 
that reenergize the afterglow emission after the jet-break 
time (at ~ 0.5 days). An alternative interpretation of 
the prominent initial brightening in terms of the emer- 
gence of the wide component in a two -co mponent out- 
flow w as given bv lBerger et al.l (|2003H) and lSheth et alJ 
(2003). They identified the early break in the opti- 
cal lightcurve with fj e t,n, the initial brightening with 
idoc,w, and a subsequent break in the radio lightcurve 
at ~ 10 days with tj e t,w- 9 To apply our two-component 
model to this source, we adopt the latter interpretation 
and consider the early afterglow lightcurve (t < 1 day) 
as having been dominated by the narrow outflow com- 
ponent. We incorporate the apparent presence of re- 
freshed shocks by taking the value of E w at tj e t,w as being 
~ 3 times larger than t he corresponding value at idcc,w 
( see IGranot et al.l l2~0031. We adopt £ , i so ,w,52/"-o ~ 30 at 
t = £jet,w on the basis of the (rather uncertain) estimates 
obtained from fitting the sizes of t he radio images of 
this afterglow at t = 24 and 83 days l|Tavlor et alJl2004t 
IGranot et alJ 2005). Using this value in equation © and 
-Eiso.w, 52/^0 ~ 10 in equation 10}, we infer 9j iV/ = 0.26 
and t? w = 8.8, respectively. Adopting p = 2.25 from the 
early-afterglow spectral fit of lWillingale et aTl 1(20041 ) . we 
extend the narrow-component's flux from tj e t, n to idcc,w 
and deduce fi = F„ w (t dcc , w ) / F^ n (t doc<w ) w 2.3. Equa- 
tions (jU, ©, and Q2J then yield Ei SOtW /E iso>B = 0.38 
and ?7 w 0j in = 0.76. We thus infer 6» j>n ' w 0.086. As a 
check on these deductions, we calculate the flux ratio 
A = ^,w(ijot,w)/^n(ijot,n) using equation We 
obtain /1 w 0.025, which agrees well with the observed 
value of ~ 0.033 given that ~ 20% of the observed flux 
at tj e t,w appears to have come from the associated super - 
nova SN 2003dh (jBerger et al.l2003btlLipkin et al.l2004]) . 
To estimate the true energy of the two outflow compo- 
nents, we assume that the measured is qtropic-equiyalent 
7-ray energy of this burst (~ 10 52 ergs: lPrice et al.ll2003t 
iHiorth et al.l 12003) has been produced by the narrow 
component with a radiative efficiency ~ 0.2, which im- 
plies £dso,n,52 ~ 5 and hence no ~ 0.19. The latter value 
is consistent with the density estimates o btained from 
spectral modeling of the afterglow (e.g., iBerger et alJ 
l2003btlWilliT 1g a,le et a,l J2004HGra,not et a,l.l2005l). Tn this 
way we deduce 2? w ,5o ~ 6.4 and -E n ,5o ~ 1-8. 

An illustrative two-component model lightcurve based 
on the above estimates is shown in Figure 01 The val- 
ues of the parameters e e and es in this fit were chosen 

9 A milder break around 0.25 days can be interp reted as cor- 
responding to the transition time t c (eq. 1121 ): see ILipkin et all 
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Fig. 3. — R-band lightcurve from a two-component model with 
parameters appropriate to GRB 030329 (Dl = 0.8 Gpc): r) n = 200, 
»?„ = 8.8, E w = 6.4 X 10 50 ergs, E n = 1.8 X 10 50 ergs, 6» jjW = 0.26, 
j n = 0.086, n = 0.19 cm -3 , e e = 0.01, eg = 0.008, and p = 2.25. 



somewhat arbitrarily to approximate the measured flux 
level (which is mostly sensitive to e e ). The value of r) B 
cannot be inferred from the available observations and 
was simply chosen to be 3> i] w . Although this simple 
model can account for the presence of the pronounced 
bump at t « 1.5 days, the model lightcurv e (which in 
general cannot rise faster than cx t ; see eqs. 

mo-Em) 

cannot readily reproduce the sharpness of the observed 
flux increase at that time. The steep rise might, how- 
ever, be explained by the refreshed-shock model (e.g., 
IGranot et al.l 12003). This possibility is consistent with 
observational evidence that the "shock refreshing" pro- 
cess in this source was already under way at t = idec,w 
ijWillingale et all 120041 iLipkin et alJ l2004[) . As noted 
above in connection with GRB 021004, the occurrence 
of such shocks around ij e t,w can be naturally expected 
in the hydromagnetic, initially neutron-rich outflow sce- 
nario. It is worth pointing out, however, that other po- 
tential problems with the two-component interpretation 
of the GRB 03032 9 data still rem ain to be addressed. In 
particular (sec Pira n et alJl2004|) . the strong radio signa- 
ture from a reverse shock that is expected in this picture 
at t w id w has not been detected. 

The afterglows of GRB 021004 and GRB 030329 were 
distinguished by the fact that their monitoring started 
early on and was conducted with particularly high time 
resolution and precision. Future observations will de- 
termine whether the departure from a smooth power- 
law behavior exhibited by the lightcurves from these two 
sources is a common feature of GRB afterglows. If this 
turns out to be the case, then the possibility that this 
behavior is associated with the decoupling of a neutron 
component in a hydromagnetically driven jet would merit 
a closer examination. 

3.2. XRF Afterglows and Source Energetics 

X-ray flash (X RF) sources (e.g., iHeTsel l2003t 
IKippen et all 12003]) arc high-energy transients that 
strongly resemble GRBs except that their peak ener- 
gies fall in the X-ray, rather than the 7-ray, spectral 
regime. One attractive interpretation of these sources 
is that they represent essentially uniform GRB jets that 



are o bserved outside the jet half-o pening angle, o b s > 9j 
(e.g.. lYamazaki et al]l2002L I2004D . In this picture, the 
larger viewing angle results in a smaller Doppler fac- 
tor and hence a lower apparent peak frequency than 
in GRBs, which correspond to 9 ^ s < 9 y The associ- 
ation with GRBs has received support from the detec- 
tion of aft erglow emission i n seve ral XRF sources. In 
particular, iSoderberg et alJ l)2004j) carried out the first 
spectroscopic observations of an XRF and form model- 
ing of the radio afterglow of XRF 020903 inferred that 
its total kinetic energy is comparable to that typically 
deduced in GRB sources. The relatively low isotropic 
equivalent energy of the prompt-emission in XRF 020903, 
Ex,iso ~ 1.1 x 10 49 ergs, is ~ 2 orders of magnitude 
smaller than the narrowly clustered values for the true 
energy out put i n gam ma-rays deduced for GR Bs, E 1 ~ 
10 51 ergs (|Frail et alJl200lt iBloom et all 12001 . which 
is itself ^1 — 3 orders of magnitude smaller than the 
isotropic equivalent gamma-ray energy output in GRBs, 
^7,iso ~ 10 52 — 10 54 ergs. This is consistent with the 
expected reduction in the measured fluence for off-axis 
observers as a result of the decrease in the Doppler fac- 
tor. 10 

In the context of a two-component outflow model with 
E w > E n , the above unified GRB/XRF picture leads 
to the identification of XRFs with GRB outflows that 
are observed at f? b s > #j. n but likely still within the 
opening half- angle of the wide component. Using our 
fiducial model parameters and assuming (9 b s = 2#j, n , 
we find that tdccwAe.n = 0.65 and 0.75 (eq. [231 ) and 
that [F v , n /F u ^](t = t e , n ) = 0.43 and 0.28 (eq. jSU) 
for E w /E n — 2 and 3, respectively. Thus, even though 
the wide component would dominate the overall after- 
glow emission, the narrow component might give rise to 
a bump in the optical lightcurve around tg n . In fact, a 
bump in the afterglow lightcur ve of XRF 030723 w as in- 
terpreted along these lines bv iHuang et al.l (|2004), who 
employed a two-component outflow model with hW « 
3#j )n and E n « 3E W . To account for the relatively late 
occurrence of this bump (between ~ 11 and ~ 14 days), 
a rather large observation angle was adopted in this fit 
(#obs = 0.37 rad w 4 0j >n ; see eq. j2Hl)- However, the re- 
brightening of XRF 030723 was foun d to be accompanie d 
by a significant spectral reddening l)Fvnbo et alJ l2004"l. 
which is not naturally explained in the two-component 
model but could possibly be associated with a supernova. 
Further tests of this aspect of the two-component model 
would therefore need to await the detection of additional 
examples of bumps in XRF afterglow lightcurves by fu- 
ture observations. 

An alternative interpretation of XRFs has been pro- 
posed in the context of the E„ > j? w two-compo nent 
collapsar outflow model (e.g.. lW. Zhang et al.ll2004j) . In 
this picture, the transient X-ray emission is attributed 
to external shocks driven by the wide (cocoon) compo- 
nent. Although the detailed implications of this proposal 
have not yet been fully worked out, this scenario could 
conceivably account for the large inferred ratio of the 
transient radiation energy and the afterglow kinetic en- 

10 The fluence dec reases as the third power of the Doppler factor 
( Gra not et al"l l20021. one power from the reduced energy of each 
photon and a power of two from the relativistic beaming of the 
photons (aberration of light) away from our line of sight. 
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ergy in XRF 020903 by associating the afterglow emission 
with the external shock of the more energetic narrow (jet- 
core) component. However, since the afterglow emission 
from the narrow component only becomes observable for 
t > te,ui it should exhibit a rapid decline with time, 
which may no t be consistent with t he data from XRF 
020903 ( where iSoderberg et al.ll2004l still detected radio 
afterglow emission after ~ 200 days). 

4. CONCLUSIONS 

In this paper we study the optical afterglow lightcurves 
produced by GRB sources that have two distinct outflow 
components. The possibility that GRB jets have a nar- 
row/fast core and a wider/slower outer component has 
been indicated by observations of both the 7-ray and 
the afterglow emission (including the afterglows of GRB 
970508, GRB 991216, and GRB 030329) and indepen- 
dently by theoretical considerations. Since we are inter- 
ested in the distinct afterglow signatures of the two com- 
ponents, we focus on models in which the Lorentz fac- 
tors of both outflow components are initially 3> 1. We 
choose the hydr omagnetically driven, initially neutron- 
rich jet model of Vlahakj^^jJJj 2003) an d the collapsar 
jet-breakout model of IW. Zhang et al.1 l)2004fl as being 
representative of scenarios in which the dual nature of 
the outflow reflects initial conditions and propagation ef- 
fects, respectively. In both of these models the charac- 
teristic initial Lorentz factor and opening half-angles arc 
rj a > 10 2 and 0j iU ~ 0.05 for the narrow component and 
rj w ~ 10, 9j iW < 'S0j lU for the wide one, and the 7-ray 
emission originates in the narrow component. They are 
distinguished, however, by the ratio of the kinetic energy 
injected into the two components: E w /E n ~ 0.1 for the 
collapsar model and > 2 for the neutron-rich hydromag- 
netic model (with E w + E n inferred to be ~ 10 51 ergs). 

Using a simple synchrotron emission model, we calcu- 
late the afterglow emission produced by the shocks that 
the two components drive into the ambient medium. We 
derive useful algebraic expressions for the component flux 
ratios at the main transition times in the light curve (in 
particular, the wide component's deceleration time tdcc,w 
and the jet-break times £j e t,w and ijet.n) for the cases 
where the observation angle £? bs satisfies £? bs < 1.5 0j, n 
and 1.5 6>j, n < 6* b s < 1.5 0j, w , respectively (where in the 
latter case t$ ttl is the relevant break time for the nar- 
row component). We study the behavior of the optical 
lightcurves for different values of E w /E n and find that, 
for the adopted characteristic parameters, the contribu- 
tion of the narrow component dominates at all times if 
this ratio is <gC 1 (as in the collapsar jet-breakout model), 
but that the contribution of the wide component be- 
comes dominant for t > £dec,w if E w > 2 E n (as in the 
neutron-rich hydromagnetic model). The emergence of 
the wide component may be related to the pronounced 
brightening detected in the lightcurves of several after- 
glows - 0.1 - 1 days after the GRB (see Fig.|3J). 

For typical parameter values idcc,w is found to be com- 
parable to ijct.n- It follows that, if E w > E n , then 



the steepening of the narrow component's lightcurve at 
t > tjet.n could be masked by the emergence (and sub- 
sequent dominance) of the wide component (see Fig. 
Under these circumstances, the only clearly discernible 
jet break in the optical lightcurve would occur at ij e t,w 
We suggest that this may have led to an overestimate of 
the emitted 7-ray energy in many GRBs because the wide 
component's opening half-angle 9^ w — rather than the 
narrow component's angle hrl — was used in convert- 
ing the measured E 7i i so into the true energy £L. This, 
in turn, would have led to an overestimate [by a fac- 
tor ~ {E n /E VJ )(9- hVJ /e- hn ) 2 = E iso , n / E iso , w ] of the ratio 
E 7 /E n that determines the kinetic-to-radiative energy 
conversion efficiency of the outflow. Factoring in this 
overestimate [which can be done when the component 
kinetic energies satisfy 1 < E w /E n < (#j,w/#j,n) 2 ] would 
alleviate the need to account for conversion efficiencies 
O(l) in internal shock models of GRBs. Dense monitor- 
ing of the afterglow lightcurve during the time interval 
(~ 0.1 — 1 days) that encompasses tj e t,n and idcc,w could 
provide a test of this suggestion. If £j t,n ^ £dcc,w then the 
lightcurve should exhibit a convex bump during this time 
interval, whereas if this inequality is reversed a concave 
flattening would be expected (see Fig.|2J. The above con- 
siderations also apply to jets observed at o bs ^ 1.5 ^j, n , 
which would be perceived as X-ray flash sources. For 
E w > E n , the afterglow emission from such sources would 
be dominated by the wide outflow component, although 
the narrow component might give rise to a bump in the 
lightcurve at t « tg <n . 

The well-monitored afterglow lightcurves of GRB 
021004 and GRB 030329 exhibited a significant early 
brightening that was followed by several less pronounced 
rebrightenings on a time scale of days. These episodes 
can be satisfactorily interpreted as refreshed shocks. We 
point out that the initially neutron-rich hydromagnetic 
jet model, in which the decoupled protons and neutrons 
give rise to a narrow/fast and wide/slow outflow compo- 
nents, respectively, could naturally account for the ap- 
pearance of refreshed shocks following (or even coincident 
with) the wide component's emergence in the afterglow 
lightcurve. Future high-quality observations should be 
able to determine whether the nonsteady behavior found 
in these two objects is a common trait of GRB sources. 
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APPENDIX 

R-BAND FLUXES FROM THE TWO OUTFLOW COMPONENTS 

For the narrow/fast jet component, typically idcc.n is smaller than all the other transition times and to,n < tjet,n> 
so there are four interesting regimes of observed frequencies i v > i^o,n 5 ^c,n(^jet,n) ^ ^ ^ ^o.n! ^m.n(^jet,n 



Two-Component Jet Models of GRB Sources 
^c,u(*jet,n); v < v m,n(tjet,a)- The corresponding fluxes are 

(^dcc,n/^c.n) ^ (^/^dcc,n) ^ 

(iAc„) 1/6 

(i/ic.„)- 1/4 

(t It )- 1/4 (t/t H 3p ~ 2)/4 
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(A3) 



(A4) 



It is worth noting at this point that the simple power-law scalings we employ in this paper may not always give an 
accurate representation of the actual behavior of a real outflow. Possible discrepancies have, in fact, been indicated 
by the hydrodynamical simul ation s reported in lGranot et alJ (j2£)01). For example, the scaling F v oc i -1 / 3 given in the 
next-to-last line of equation (|A4|) appears to differ from the behavior exhibited in Figure 2 of that reference, where 
the flux at low frequencies is shown to continue rising well past tj e t. However, the expressions that are most relevant 
to the behavior of the optical afterglow seem to be consistent with the results of the numerical simulations. 

For the wide/slow jet component, the cooling frequency v c is always larger than the characteristic frequency v m 
for typical parameters, so it is always slow cooling. This can be verified in the following analysis. Since the ratio 
Vcjvm decreases with time before idee and increases with time after idee, it reaches its minimum value at idee, and the 
condition for always being in slow cooling is {v c /v m )\t dac > 1- Combining equations Q and (J2J), we get 



{v c lv m )\ Ua 



oo„-2 p -2/3 -2 -2 
O.oy - c 'i S o,w,52 t B,-l fc e,- 



(A5) 



For the typical parameter ranges described at the end of this Appendix, {v c /v m )\t dcc ~ 5 x 1CP 5 — 2 x 10 9 , which shows 
that (t / c /f m )|t doc w > 1 for most parameter values. This inequality is violated only when all the relevant parameters are 
close to their maximum values, which is not a typical situation. If, in addition, f c ,w(ijet,w) > ^m,w(idcc,w) is also satisfied 
(which is again true for most parameter values), then there are five interesting frequency regimes: v > f c . w (idcc,w), 

^c,w(^jet,w) *\ ^ ^ ^cw (idcc,w) 5 ^ / m,w(idec,w) *\ ^ *C ^c,w(^jet,w)> ^m,w(ijct,w) ^ ^ ^ ^m,w(^dec,w)j and V <C ^m,w(ijct,w)- 

The corresponding fluxes are 
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If, on the other hand, ^ c ,w(^jet,w) < ^m,w(^dec,w)) then there is an additional possible shape for the lightcurve when 
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In this work we are interested in the R-band spectral regime, so we compare the characteristic frequencies at the 
main spectral transition times with a typical optical frequency i/r. In particular, at £dec,w> 



^CjW^dcc^w) 
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We also have 



* c (i jet ) = 1.6 x 10 14 B^e-^n- 5 / 6 ^ 3 Hz 
i^(t jet ) = 1.4 x 10 13 g 2 e 2 e 
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(A12) 
(A13) 
(A14) 

(A15) 

(A16) 
(A17) 



We adopt as typical parameter ranges p ~ 1.5 — 3, no = 0.3 — 30, e e = 0.005 — 0.3, e_e = 0.001 — 0.1, 9j yU ~ 
0.05 - 0.1, jlW - 0.1 - 0.3, E n ~ E w ~ 10 50 - 10 51 ergs l|Panaitescu fc Kumarl 12002]) . 7? n ~ 10 2 , and r] w ~ 10.' For 
this set of parameters we find that the following inequalities are always obeyed: ^m,n(idec,n) > > ^m,n(^jot,Ti), 
z/r > ^ m ,n(idcc,w), and z/r > f m . w (fdcc,w) • The observation frequency z/r can be larger or smaller than v c for both 
components at tj e t,n and tdec.w However, for the parameter combinations employed in the plots shown in this paper, 
^c,n(tjet,n) > vr > ^m,n(*jet,n) and ^ c ,w(idcc,w) > vr > ^m,w(idcc,w)- Since v m and i/ m /i/ c decrease with time (see 
eqs. and [S]), the inequalities v c > z^ m and z/r > v m continue to apply after the specified times. For the narrow 
component v c>n > z/r > z/ ni!n also continues to apply since (under the assumptions underlying eq. [S]) z/ Cin remains 
constant for t > tj C t,n- However, z/ CjW decreases with time after idcc.w and could potentially fall below z/r before its 
value becomes frozen at t\_ w . The foregoing arguments imply that the frequency regimes corresponding to equations 
<|A4I) . (|X9|) . l|A10jl . and IjAlljl would typically not be relevant. 



Two-Component Jet Models of GRB Sources 



13 



REFERENCES 



Beloborodov, A. M. 2000, ApJ, 539, L25 

Berger, E., Kulkarni, S., & Frail, D. A. 2003a, ApJ, 590, 379 
Berger, E., et al. 2003b, Nature, 426, 154 
Blandford, R. D., & McKee, C. 1976, Phys. Fluids, 19, 1130 
Bloom, J. S., Frail, D. A., & Kulkarni, S. R. 2003, ApJ, 594, 674 
Cannizzo, J. K., Gehrels, N., & Vishniac, E. T. 2004, ApJ, 601, 
380 

Fox, D. W., et al. 2003, Nature, 422, 284 
Frail, D. A., et al. 2000, ApJ, 538, L129 
Frail, D. A., et al. 2001, ApJ, 562, L55 
Fynbo, J. P. U., et al. 2004, ApJ, 609, 962 

Granot, J., Miller, M., Piran, T., Suen, W. M., & Hughes, P. A. 

2001, in Gamma-Ray Bursts in the Afterglow Era, ed. E. Costa, 
F. Frontera, & J. Hjorth (Berlin: Springer), 312 

Granot, J., Nakar, E., & Piran, T. 2003, Nature, 426, 138 
Granot, J., Panaitescu, A., Kumar, P., & Woosley, S. E. 2002, ApJ, 
570, L61 

Granot, J., Ramirez-Ruiz, E., & Loeb, A. 2005, ApJ, 618, 413 
Heise, J. 2003, in AIP Conf. Proc. 662, Gamma-Ray Burst and 

Afterglow Astronomy 2001, ed. G. R. Ricker & R. K. Vanderspek 

(Melville: AIP), 229 
Hjorth, J., et al. 2003, Nature, 423, 847 

Huang, Y. F., Wu, X. F., Dai, Z. G., Ma, H. T., & Lu, T. 2004, 
ApJ, 605, 300 

Khokhlov, A. M., Hofiich, P. A., Oran, E. S., Wheeler, J. C, Wang, 

L., & Chtchelkanova, Y. 1999, ApJ, 524, L107 
Kippen, M., Woods, P. M., Heise, J., In't Zand, J. J. M., Briggs, M. 

S., & Preece, R. D. 2003, in AIP Conf. Proc. 662, Gamma-Ray 

Burst and Afterglow Astronomy 2001, ed. G. R. Ricker Sz R. K. 

Vanderspek (Melville: AIP), 244 
Kobayashi, S., & Sari, R. 2001, ApJ, 551, 934 
Kobayashi, S., k, Zhang, B. 2003, ApJ, 582, L75 
Konigl, A. 2004, in AIP Conf. Proc. 727, Gamma-Ray Bursts: 30 

Years of Discovery, ed. E. E. Fenimore & M. Galassi (Melville: 

AIP), 257 

Kumar, P., & Granot, J. 2003, ApJ, 591, 1075 
Kumar, P., & Piran, T. 2000a, ApJ, 532, 286 
Kumar, P., & Piran, T. 2000b, ApJ, 535, 152 

Lazzati, D., Rossi, E., Covino, S., Ghisellini, G., & Malesani, D. 

2002, A&A, 396, L5 

Levinson, A., & Eichler, D. 1993, ApJ, 418, 386 
Liang, E. W., & Dai, Z. G. 2004, ApJ, 608, L9 
Lipkin, Y. M., et al. 2004, ApJ, 606, 381 
Meszaros, P. 2002, ARAA, 40, 137 
Nakar, E., & Oren, Y. 2004, ApJ, 602, L97 
Nakar, E., & Piran, T. 2002, ApJ, 572, L139 



Nakar, E., & Piran, T. 2003, ApJ, 598, 400 

Nakar, E., Piran, T., & Granot, J. 2002, ApJ, 579, 699 

Nakar, E., Piran, T., & Granot, J. 2003, New A, 8, 495 

Panaitescu, A., & Kumar, P. 2002, ApJ, 571, 779 

Panaitescu, A., Meszaros, P., & Rees, M. J. 1998, ApJ, 503, 314 

Pedersen, H., et al. 1998, ApJ, 496, 311 

Piran, T. 1999, Phys. Rep. 314, 575 

Piran, T., Nakar, E., & Granot, J. 2004, in AIP Conf. Proc. 727, 
Gamma-Ray Bursts: 30 Years of Discovery, ed. E. E. Fenimore 
& M. Galassi (Melville: AIP), 181 

Price, P.A., et al. 2003, Nature, 423, 824 

Ramirez-Ruiz, E., Celotti, A., & Rees, M. J. 2002, MNRAS, 337, 
1349 

Rees, M. J., & Meszaros, P. 1994, ApJ, 430, L93 
Rhoads, J. 1999, ApJ, 525, 737 

Rossi, E., Lazzati, D., & Rees, M. J. 2002, MNRAS, 332, 945 

Sari, R. 1997, ApJ, 489, L37 

Sari, R., & Meszaros, P. 2000, ApJ, 535, L33 

Sari, R., & Piran, T. 1995, ApJ, 455, L143 

Sari, R., & Piran, T. 1999a, ApJ, 520, 641 

Sari, R., & Piran, T. 1999b, A&A Supl. Ser., 138, 537 

Sari, R., Piran, T., & Halpern J. P., 1999, ApJ, 519, L17 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Sheth, K., et al. 2003, ApJ, 595, L33 

Soderberg, A. M., et al. 2004, ApJ, 606, 994 

Starling, R. L. C, Wijers, R. A. M. J., Hughes, M. A., Tanvir, N. 

R., V reeswijk, P. M., RoL E., & Salamanca, I. 2005, MNRAS, in 

press lastro-ph/0501120l 
Taylor, G.B., Frail, D.A., Berger, E., & Kulkarni, S.R. 2004, ApJ, 

609, LI 

Uemura, M., Kato, T., Ishioka, R., & Yamaoka, H. 2003, PASJ, 
55L, 31 

van Putten, M. H. P. M., & Levinson, A. 2003, ApJ, 584, 937 
Vlahakis, N., Peng, F., & Konigl, A. 2003, ApJ, 594, L23 
Willingale, R., Osborne, J. P., O'Brien, P. T., Ward, M. J., Levan, 

A., & Page, K. L. 2004, MNRAS, 349, 31 
Yamazaki, R., Ioka, K., & Nakamura, T. 2002, ApJ, 571, L31 
Yamazaki, R., Ioka, K., & Nakamura, T. 2004, ApJ, 606, L33 
Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 

597, 459 

Zhang, B., Dai, X., Lloyd-Ronning, N. M., & Meszaros, P. 2004, 

ApJ, 601, L119 
Zhang, B., & Meszaros, P. 2002, ApJ, 571, 876 
Zhang, W., Woosley, S. E. & Heger, A. 2004, ApJ, 608, 365 



